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Abstract. Fractional diffusion equations replace the integer-order derivatives in 
space and time by their fractional-order analogues. They are used in physics to 
model anomalous diffusion. This paper develops strong solutions of space-time 
fractional diffusion equations on bounded domains, as well as probabilistic repre- 
sentations of these solutions, which are useful for particle tracking codes. 



1. Introduction 

The traditional diffusion equation dtu = Au describes a cloud of spreading par- 
ticles at the macroscopic level. The point source solution is a Gaussian probability 
density that predicts the relative particle concentration. Brownian motion provides a 
microscopic picture, describing the paths of individual particles. A Brownian motion, 
killed or stopped upon leaving a domain, can be used to solve Dirichlet boundary 
value problems for the heat equation, as well as some elliptic equations [U |TT]. The 
space-time fractional diffusion equation d^u = A"/^m with < /3 < 1 and < a < 2 
models anomalous diffusion [19]. The fractional derivative in time can be used to 
describe particle sticking and trapping phenomena. The fractional space derivative 
models long particle jumps. The combined effect produces a concentration profile 
with a sharper peak, and heavier tails. This paper studies strong solutions, and prob- 
abilistic representations of solutions, for the space-time diffusion equation on bounded 
domains. Our main result is Theorem 15.11 Strong solutions are obtained by separa- 
tion of variables, combining the Mittag-Leffler solution to the time-fractional problem 
with an eigenfunction expansion of the fractional Laplacian on bounded domains. The 
probabilistic representation of solutions involves an inverse stable subordinator time 
change, resulting in a non-Markovian process. Fractional diffusion equations are be- 
coming popular in many areas of application [151 [23]. these applications, it is often 
important to consider boundary value problems. Hence it is useful to develop solu- 
tions for space-time fractional diffusion equations on bounded domains with Dirichlet 
boundary conditions. 
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2. Random walks and stable processes 



A random walk St = Yi + ■ ■ ■ + Y[t], a. sum of independent and identically distributed 
M'^-valued random vectors, is commonly used to model diffusion in statistical physics. 
Here [t] denotes the largest integer not exceeding t, and Sn represents the location of 
a random particle at time n. Suppose the distribution of Y is spherically symmetric. 
If (T^ := E[|yip] is finite and ]E[Yi] = 0, Donsker's invariance principle implies that 
as A — oo, the random process {X~^^'^Sxt, t > 0} converges weakly in the Skorohod 
space to a Brownian motion {Bt, t > 0} with IE[-Bi] = cr^- If the step random variable 
Yi is spherically symmetric, and P(|Yi| > x) ~ as x — > oo for some < a < 2 

and C > 0, then E is infinite, and the extended central limit theorem tells us 

that {\~^/''Sxt,t > 0} converges weakly to a rotationally symmetric a-stable Levy 
motion {At,t > 0} with 

E[e'^-^*] = e-^°l«l"* for every ^ G M'^ and t > 0, 

where the constant Co depends only on C and the dimension d, see [18]. A simple 
rescaling in space yields a standard stable process with Cq = 1. Since {X^''°'At,t > 0} 
has the same distribution as {Axt,t > 0}, stable Levy motion represents a model 
for anomalous super- diffusion, where particles spread faster than a Brownian motion 

If we impose a random waiting time T„ before the nth random walk jump, then 
the position of the particle at time T„ = Ji + ■ ■ ■ + J„ is given by Sn- The number 
of jumps by time t > is iVf = max{n : T„ < t}, so the position of the particle at 
time t > is SNt, a subordinated process. If P(J„, > t) ~ Ct~^ as t — )■ oo for some 
< /3 < 1, then the scaling limit of c~^^^T[ct] ^ as c — oo is a strictly increasing 
stable Levy motion with index /3, sometimes called a stable subordinator. The jump 
times T„ and the number of jumps Nt are inverses: {Nt > n} = {T„ < t}. [20} 
Theorem 3.2] shows that {c'^^Nctjt > 0} converges weakly to the process {Et,t > 0}, 
where Et = mi{x : Z^j. > t}. In other words, the scaling limits are also inverses: 
{Et < x} = {Zr, > t}. Now Net ~ c'^Et, and [201 Theorem 4.2] shows that the scaling 
limit of the particle location {c'^/^S'atj^^j, t > 0} is {AEt,t > 0}, a symmetric stable 
Levy motion time-changed by an inverse stable subordinator. 

The random variable Zt has a smooth density. For properly scaled waiting times, 
the density of the standard stable subordinator Zt has Laplace transform E[e~''^'] = 
Q-tv'^ for any ri,t > 0, and Zt is identically distributed with t^^l^Zi. Writing g^i^u) for 
the density of Zi, it follows that Z<, has density s~^^^gi^{s~^^^u) for any s > 0. Using 
the inverse relation F{Et < s) = F{Zs > t) and taking derivatives, it follows that Et 
has the density 

(2.1) Ms) = ^nZs >t)= 0~h~'~'/^gp{ts~'/^). 

as 

For more details, see [191 120] • 
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3. Fractional calculus 



The Caputo fractional derivative of order < /3 < 1, defined by 

d^f{t) 1 f'dfir) dr 



was invented to properly handle initial values [HI [I3]. Its Laplace transform (LT) 
s^f{s) — s^~^f{0) incorporates the initial value in the same way as the first derivative. 
Here f{s) = e~^^f{t)dt is the usual Laplace transform. The Caputo derivative has 
been widely used to solve ordinary differential equations that involve a fractional time 
derivative [13 123] . In particular, it is well known that the Caputo derivative has a 
continuous spectrum, with eigenfunctions given in terms of the Mittag-Leffler function 

OO JU 



In fact, f{t) = Eii{—\t^) solves the eigenvalue equation 

^'f^'^- Xf(t) 

for any A > 0. This is easy to check, differentiating term-by-term and using the fact 
that fP has Caputo derivative tP-^T{p + 1) /T{p 1 - /3) for p > and < /3 < 1. 
For < a < 2, the fractional Laplacian A°/^/ is defined for 

/GDom(A"/2) = //eL2(M'^;dx): / \^nm\'d^ < oo 

as the function with Fourier transform 

(3.2) A^(0 = -ier7(0- 

For suitable test functions (for example, functions with bounded second deriva- 
tives), the fractional Laplacian can be defined pointwise: 

(3.3) A"/V(x) = / {fix + y)- fix) - Vfix) ■ y^\y\<^}) ^ dy, 

Jym'^ \y\ 

where Q,a > is a specific constant that depends on d and a so that 

f l-cos?/i 



Remark 3.1. (i) It can be verified using Fourier transforms that, for / E Dom(A"/^), if 
the right hand side of (13.31) is well-defined for a.e. x G M"', then the Fourier transform of 
the right-hand side of (13.31) equals — |'C|"/(0 (cf- [HI Theorem 7.3.16]). Conversely, it 
can also be verified that if / G L^(]R'^; dx) is a function such that the right hand side of 
(13.31) is well-defined for a.e. x eW^ and is L^(]R'^; (ix)-integrable, then / G Dom(A°/^) 
and (13.31) holds. 
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(ii) Using a Taylor series expansion in f l3.3p . it is easy to see that A"/^/(a;o) exists 
and is finite at a point xq G M'^ if / is bounded on M.'^ and / is at tlie point xq- 
Hence, if / is bounded and continuous on M'^ and / is in an open set D, then A"/^/ 
exists pointwise and is continuous in D. Moreover, if / is a function on [0, oo) 
with |/'(^)| ^ cf^'^ for some 7 > 0, then by (13.1 1) , the Caputo fractional derivative 

f{t)/dt^ of / exists for every t > and the derivative is continuous in t > 0. □ 

For < a < 2, let X be the Levy process on such that 

This Levy process X is called a standard (rotationally) symmetric a-stable process 
on M'^. When a = 2, it is Brownian motion running at double speed. 

Denote the transition semigroup of X by {P^, t > 0}. Using the fact that Xf =^ Xq 
as t ^ 0+, it is not hard to show (e.g., see [H Theorem 13.4.2]) that {Pt,t > 0} 
is a symmetric strongly continuous semigroup on the Banach space L^{M.'^; dx). Let 
(J", £) be the Dirichlet form of X on ^^(IR'^; dx). That is, 

(3.4) J" = \u e L^iR'^^dx) : suphu- PtU,u)L2(Rd.dx) < oo\ , 

L t>o t ) 

(3.5) £{u,v) = - PtU,v)L2(^d-dx) ioi u,v e 

t-s>0 t 

It is known that, for example, via Fourier transforms [1 



^ = := L G L\R''; dx) : [ ^^L^M}ldxdy < 00 

I jR'^xR'i \x ~ y\ " 

£{u,v) = ^ / {u{x)-u{y)){v{x)-v{y))^^^^_ 

Let (Dom(£), £) be the L^-generator of the Dirichlet form {£, J^); that is, / G Dom(£) 
if and only if / G W°'^'^''^{M.'^) and there is some u G L^(M'^; c/x) so that 

S{f, g) = -(n, g) for every ^ G PF"/2,2^^d^). 

in this case, we denote this u by Cf . It is known (cf. [H]) that £ is also the semigroup 
generator of {Pt^t > 0} on the space L'^{M.'^; dx). Using the Fourier transform, one 
can conclude (cf. [H]) that / G Dom(£) if and only if f^^ |C|"|/(OP'^^ < ^5 
^fiO = every / G Dom(£). Hence the L^-generator of X is the 

fractional Laplacian A"/^. 

It follows directly from Dirichlet form theory (cf. [H]) that, for / G L^(]R'^) and 
t > 0, Ptf e T = iy"/2'2(M'i), and w(t,x) := E^[f{Xt)] is a weak solution to the 
following parabolic equation: 

(3.6) ^v{t, x) = A^/Mt, x); v{0, x) = f{x). 
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That is, the function x t— )■ v{x, t) belongs to the domain of the generator C = A"/^ 
for every t > 0, and equation fl3.6p holds in the space L^(]R'^; dx). Here the fractional 
Laplacian and the first time derivative in (13.61) are defined in terms of the Banach 
space norm. For example, the time derivative is the limit of a difference quotient 
that converges in the sense, so it need not exist point-wise. The classical diffusion 
equation models the evolution of particles away from their starting point, due to 
molecular collisions. The space-fractional diffusion equation f l3.6p models particle 
motions in a heterogeneous environment, where the probability of long particle jumps 
follows a power law [T7] . 

For < a < 2, the symmetric a-stable process X can be obtained from Brownian 
motion on M through subordination in the sense of Bochner [8]. Let {5,P^,x G M."^} 
be Brownian motion on M*^ with Wx{Bq = x) = 1 and Eo[-Bf-Bf] = 2tl, where ' denotes 
the transpose, and / is the d x d identity matrix. For < a < 2, let Zt be a standard 
stable subordinator with Zq = 0, whose Laplace transform is E[e~'^^*] = e"*""''^ for 
every s,t > 0. Then it is easy to verify, using Fourier transforms and a simple 
conditioning argument, that Bzt is a symmetric a-stable Levy process starting from 
the origin that has the same distribution as X, with Xq = 0. The process X has a 
jointly continuous transition density function p{t,x,y) = pt{x — y) with respect to 
the Lebesgue measure in M*^. That is, 

P,(XtGA)= / p{t,x,y)dy. 

J A 

Using the self-similarity of the stable process and its relation with Brownian motion 
through subordination, it is not hard to show that for a G (0, 2) we have 

(3.7) pt{x) = r^/>i(t-i/"x) < t"'^/>i(0) =: r'^/'^Md^a, t > 0,x G M"^. 

Another kind of time change relates to particle waiting times. Suppose {Tt,t > 0} 
is a uniformly bounded strongly continuous semigroup on a Banach space E, with 
infinitesimal generator (^, Dom(^)). It is known that v{t) = Ttf solves the Cauchy 
problem dv/dt = Av with t'(O) = / for any / G Dom(^) (see [2]). Let Z be a standard 
/9-stable subordinator independent of X, and recall that Ef = mi{s > : Zg > t} is 
its inverse process. If Q/siu) is the density of Zi, then [31 Theorem 3.1] shows that 
another subordinated semigroup 

POO 

(3.8) Rtf = / gp{u)T(^t/u)^f du 

Jo 

yields solutions to the time-fractional Cauchy problem: w{t) = Rtf solves 

Q^wit) = Aw; w{0) = f 

on the Banach space E for any / G Dom(^). Applying this to the transition semi- 
group {Pt,t > 0} of the symmetric a-stable process X on the space L^(M'^; dx), one 
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sees that the process Yt = can be used to solve the space-time diffusion equation 
on M'^; that is, w{t,x) = Ej.[/(1^)] is a weak solution for 

(3.9) at^"^^^' = ^"^'^(^' ^(^' 0) = /(x). 

That is, the function x H- t) belongs to the domain of the generator C = A"/^ 
for every t > 0, and equation fl3.9p holds in the Banach space L^(]R'^; dx). 

4. ElGENFUNCTION EXPANSION FOR BOUNDED DOMAINS 

Let D be a bounded open subset of M"^. Recall that X is a standard spherically 
symmetric stable process on M"^, and define the first exit time 

Td = inf{t > : Xt ^ D}. 

Let denote the process X killed upon leaving D; that is, X[^ = Xt for t < td and 
Xj^ = d ioi t > Td. Here 9 is a cemetery point added to D. Throughout this paper, 
we use the convention that any real-valued function / can be extended by taking 
f{d) = 0. The subprocess X^ has a jointly continuous transition density function 
Poit, X, y) with respect to the Lebesgue measure on D. In fact, by the strong Markov 
property of X, one has for t > and x,y E D, 

(4.1) pDit,x,y) =p{t,x,y) -E-,[p{t - m, X^^,y);TD < t] <p{t,x,y). 
Denote by {P^,t > 0} the transition semigroup of X^, that is 

P,^/(x)=E.[/(Xf)]= / pnit,x,y)fiy)dy. 

JD 

The proof of the following facts can be found in [H]: The operators {-Pj^, t > 0} form 
a symmetric strongly continuous contraction semigroup in L'^{D;dx). Let {£^,J^^) 
denote the Dirichlet form of X^, defined by ([S3])-([i5D but with {Pf,t > 0} in 
place of {Pt,t > 0}. Then J^^ is the -y/^-completion of the space C^{D) of 
smooth functions with compact support in D, denoted by Wq''^{D) in literature. Here 
Si{u,u) = S{u,u) + J^au{xYdx. Moreover, S^{u,v) = S{u,v) for M,f G 
Let Cd be the L^-infinitesimal generator of {S^,J^^); that is, its domain Dom(££)) 
consists all / e PF^''^'^(-D) such that 

^^(/> g) = -{u, g)mD;dx) for every g G W^^'^''^{D); 

for some u G L'^{D; dx); in this case, we denote this u by Cof- It is well-known (cf. 
[Hj) that Cd is the L^-generator of the strongly continuous semigroup {Pf,t > 0} 
in L'^{D;dx). For every / G L'^{D;dx) and t > 0, P^^ f G BomiCo) C Wo^^'^{D). 
Moreover u{t,x) := Pff{x) is the unique weak solution to 

du 

^ = ^^" 

with initial condition u{0,x) = f{x) on the Banach space L'^{D] dx). 
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Note that the transition kernel p£){t,x,y) is symmetric and strictly positive with 

(4.2) PD{t,x,y) <pit,x,y) <t~''^''Md,a, x,yeD,t>0 

in view of (13. 7p . In particular, one has sup^.g^ J^p(t,x,y)^dy < oo for every t > 0. 
Thus for each t > 0, is a Hilbert-Schmidt operator in L'^{D; dx) so it is compact. 
Therefore there is a sequence of positive numbers < Ai < A2 < A3 < ■ ■ ■ and an 
orthonormal basis {ipn-, n > 1} of L'^{D; dx) so that Pfipn = ^"^"Vn in L'^{D; dx) for 
every n>l and t > 0. Since for every / G L'^{D; dx), f\x) = J2'^=i {fy'4'n)'ipn{x), we 
have 

00 00 

(4.3) P,^/(x) = ^ {f,i:n)PFM^) = J2e-'-\f,i:n)M^)- 

n=l n=l 

That is, the transition density 

00 

(4.4) poit, x,y) = J2 e-^"'Mx)My)- 

n=l 

It follows from [3, Theorem 2.3] that for any bounded open subset D of R"', one has 

(4.5) cin"/'^ < A„ < can'*/'^ for every n>l. 
Using the spectral representation, one has 

(4.6) BomiCo) = |/ G L\D) : \\Cnf\\h^n) =Y.^l{f,^nf < 00 
and 

00 

^ofix) = - ^ \n{f, V'n)V'n(a;) for / G Dom(£z)). 

n=l 

For any real valued function : R — )■ M, one can also define the operator 4>{Cd) as 
follows: 

Dom(0(£^)) = |/GL2(Z^;cix):^</.(A„)2(/,V^„)2<oo 

00 

<P{CD)f = Y.<P{Xn){f,4^n)A 



n=l 



n=l 



In next section, the operator £^ defined using (j)(t) = will be utilized. 

The generator £d is also called the fractional Laplacian on D with zero exterior 
condition, denoted as A^/^l/j. We now record a lemma that gives an explicit expres- 
sion oi Cd. 



Lemma 4.1. For f G , if 

(4.7) := lim / {f{y) - f{x)) '%^J y 

exists and the convergence is uniformly on each compact subsets of D and cf) G 
L^{D;dx), then f G Dom(££i) and (p = C^f . In particular, if f is a hounded function 
m n C^ip), then f G Dom(£D) and 

^Df{x) = lim/ {f{y)-f{x)) . ciy 



:|?/-a,'|>e} 

{fix + y)- fix) - V/(x) ■ yl{\y\<i}) j^^^y- 

Proof Suppose that / G J^^ and that (f> defined by (14. 7p converges locally uniformly 
in D and is in L^{D;dx). Then for every g G C^iD), by the expression of £^{f,g) 
and the symmetry, 

S^'iLg) = \l ifix)-fiy))i9ix)-giy)) '%^J xdy 
jR-ixw \x — y\ 

= J lim/ ifiy) - fix))igiy) - gix))- — ^"^'^ dxdy 

^ J{(x,s/)GK'*xRrf:|a;-y|>e} W — X\ 

= -lim / (/ ifiy) -fix))- — ^"^'^ dy) g{x)dx 



(j)ix)gix)dx. 

Since C^iD) is f^f'-dense in W(^''^'^(_D), this implies that / G Dom(££)) and £d/ = 
on D. 

Assume now that / is a bounded function in J^^CiC^iD). Using a Taylor expansion, 
one easily sees that 

/ \fix + y) - fix) -Vfix) ■yl{\y\<i}\-^^^dy < oo for every x G D 
Jyewi \y\ 

and the integral is a continuous function on D. Set 

i^ix) = {fix + y)- fix)-Vfix)-yl[\y\<i})-^^^dy for x G £>. 

JyeR'' \y\ 

For any compact subset K oi D, let 

Ke := {2 G M"^ : there is some x G K so that |z — x| < e}. 

Defining 

IIDVIIoo = max 

l<i,j<d 



dxidxj 



8 



we have 



hm sup 



hm sup 



{f{y)-f{x))- 



{yeR'i:\y-x\>e} 



X 



d+a 



dy — ip{x] 



■.\y-x\<e} 



(/(a: + Z/)-/(x)-V/(x)-yl{|,|<i}) 



Cd,a 
7, 1 d+a 



dy 



< hm 

e-5>0 



sup \\D'f 

\y-x\<e} z(^K^ 



2X11 I |2 

oo \y\ 



Cd,i 



\y 



\d+a 



dy 



0. 



By what we have shown in the first part, this imphes that / G Dom(££)) with 
^of = "ip, which completes the proof of the lemma. □ 

The main purpose of this paper is to investigate the existence of strong solution to 
the following equation: 

90 



(4.8) 



u{t,x) = A"^^u{t,x)] xeD,t>0 



u{t,x) = 0, X e t > 0, 

u{0, x) = f{x), X E D. 

Let Coo(-D) denote the Banach space of bounded continuous functions on that 
vanish off D, with the sup norm. 

Definition 4.2. (i) Suppose that / G L'^{D]dx). A function u{t,x) is said to be a 
weak solution to (14.81) if u{t,-) G Wq''^{D) for every t > 0, \imtiou{t,x) = f{x) a.e. 
in D, and /dt^u{t, x) = A°'^'^u{t, x) in the distributional sense; that is, for every 
G Ci([0,oo) and 0G C2(Z}), 

9V(t) 



u{t, x) 



dt (j){x)dx 



8''{u{t,-),(t))i){t)dt. 



'0 / ^0 

(ii) Suppose that / G C{D). A function u(t, x) is said to be a strong solution (14. 8 p 
if for every t > 0, u{t,-) G Coo(-D), A"/^M(t, ■)(x) exists pointwise for every x E D 
in the sense of (13. 3p . the Caputo fractional derivative d^u{t,x)/dt^ exists pointwise 
for every t > and x G -D, /dt^u{t^x) = A"/^-u(t,x) pointwise in (0, oo) x D, and 
lim^o w(t, a;) = f{x) for every x E D. 

A boundary point x of an open set D is said to be regular for D iiFr^[TD{X) = 0] = 1. 
A sufficient condition for xq G dD to be regular for D is that D satisfies an exterior 
cone condition at xq, that is, there exists a finite right circular open cone V = V^q 
with vertex xq such that V^q C D'^ (cf. [lOl Theorem 2.2]). An open set D is said to be 
regular if every boundary point of D is regular for D. Assume now that D is a regular 
open set. Then [TOl Theorem 2.3] shows that {PP it > 0} is a strongly continuous 
(Feller) semigroup on the Banach space Coo{D) of bounded continuous functions on 
that vanish off -D, with the sup norm. Moreover, {PP it > 0} has the same 
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set of eigenvalues and eigenf unctions on C^{D) as on L'^{D]dx): Pfipn = e~^"Vn 
in Coq{D) (see [TUl Theorem 3.3]). In particular, every eigenfunction ipn oi the L^- 
generator Co is a bounded continuous function on D that vanishes continuously on 
the boundary dD. 

5. Space-time fractional diffusion in bounded domains 

In this section, we prove strong solutions to space-time fractional diffusion equa- 
tions on bounded domains in M'^. We give an explicit solution formula, based on the 
solution of the corresponding Cauchy problem. The basic argument uses an eigen- 
function expansion of the fractional Laplacian on D, and separation of variables. The 
probabilistic representation of the solution is constructed from a killed stable pro- 
cesses, whose index corresponds to the fractional Laplacian, modified by an inverse 
stable time change, whose index equals the order of the fractional time derivative. 

Recall that X is a rotationally symmetric a-stable process in M.'^ and {Et,t > 0} 
is the inverse of a standard stable subordinator of index /3 G (0, 1), independent of 
X. In the following proof, we denote by c, Ci, C2, . . . a constant that may change from 
line to line. 

Theorem 5.1. Let D be a regular open subset o/M'*. Suppose f G Dom(£|,) for 
some k > -1 + {3d + 4) / {2a) . Then 

u{t,x) = E,[/(X|J] G a([0,oo) X M'^) n Ci'2((0,oo) x D) 

andu{t,x) is a strong solution to the space-time fractional diffusion equation fl4.8p . 

Proof. First we will prove that / G Coo{D). Let < Ai < A2 < A3 < ■ ■ ■ be the 
eigenvalues of Cd and {ipnin > 1} be the corresponding eigenfunctions, which form 
an orthonormal basis for L^{D\ dx). Note that, since D is a regular open set, we have 
from the last section that ipn ^ C*oo(-D) for each n > 1. Since / G Dom(£^) for some 
> -1 + (3(i + 4) /(2a), using (USD it follows that 

00 

(5.1) M:=5^A^'=(/,^„)2<oo, 

n=l 

and so |(/,^/^„,)| < a/MA~^ From (H^]) and (i3D we get 

00 

e-""*|^A„(x)|2 < Y,^-^''\M^)?=PD{t,x,x) < M,,,r'^/" 

k=l 

and hence, taking square roots of both sides, we get 

|^„(x)| <e^"*/yM,,,t-'^/" 

Taking t = 1/A„ gives us 

(5.2) \'^n{x)\ < cA^/(2") for every xeD 
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for some c > 0. Since k > —1 + {3d + 4)/ (2a), (15 ■2p together with (14 .Sp imphes that 

CO oo oo 

n=l n=l n=l 

Hence /(x) = X]n=i {f^i^n)'ipn converges uniformly on D, and so / G Coo(-D)- 

Recall that PP f{x) = E,x[f{XP)] is the unique weak solution in Wq^'^''^{D) of the 
equation 

(5.3) ' 



dt 



v{t, x) = A°/\(t, x) with w(0, x) = f{x) 



on the Banach space L'^{M.'^; dx) (cf. (see [H]). The semigroup Pf has density 
function pD(t,x,y) given by (14. ip . Note that p(t,x,y) is smooth in x. By a proof 
similar to [H Proposition 3.3], we have for every j > 1 and 1 < i < d that 

9^' 



(5.4) 



dx: 



< c t-^-^+J')/" A 



t 

\x - 



< cit-'/yit,x,y). 



In view of the symmetry p(t,x,y) = p(t,y,x) and pD(t,x,y) = pDit,y,x), we have 
from (14. ip and (15. 4p that Pff{x) = J^p£)it,x,y)f{y)dy is smooth in x G -D. More- 
over, for every compact subset K of D and T > 0, there is a constant C2 = C2((i, a, T) 
such that, for x G -fC and t G (0, T], 



(5.5) 



dx] 



■PD{t,x,y) < Cat ^/Xt,x,y). 



The Chapman- Kolmogorov equation implies 



p{t,x,yfdy 



pit, X, ?/)p(t, x) dy = p{2t, X, x). 



It then follows using (14. 2p . (15. 5p . and the Cauchy-Schwarz inequality that 



(5.6) 



\v^prf{x)\<csr^/^{2ty 



■d/{2a) 



Consequently, each eigenfunction ipnix) = e^"'^Pftpn{x) is smooth inside D with 



V^'V'„(x) 



-{d+2j)/{2a) X^t 



for X G -ftT and t G (0,T]. Taking t = 1/A„ yields 



(5.7) 



VV„,(x) 



< C3A(f+2,)/(2a) 



for X E K. 



In view of (14. 3p . Pf f{x) is also differentiable in t > 0. (The eigenfunction expansion 
(14. 3 p together with (15. 7p gives another proof that P^^ f is C°° in x G -D.) Hence in 
view of Remark [3. 1[ v{t,x) = P^ f{x) is a classical solution for dv/dt = Cdv in D. 



Now define 



u{t,x) = E,[/(X^J] = E,[/(X^J;E, < m]. 
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Since generates a strongly continuous (Feller) semigroup on C^{D), Pff{x) is a 
bounded continuous function on [0, oo) x M.'^ that vanishes on [0, oo) x D'^, and hence 
so is u, in view of (13. 8p . By [3l Theorem 3.1] (and [201 Theorem 4.2]), u{t,x) is a 
weak solution for the parabolic equation (14. 8 p on L'^(M.'^, dx). Then, to show that u is 
a classical solution, by Remark |3.H it suffices to show that u{t, ■) G C^{D) for each 
t > 0, and that the Caputo derivative of t i— )■ u(t,x) exists for each x, and is jointly 
continuous in (t,x). 

Bingham [B] showed that the inverse stable law Et with density ft{s) given by (12.11) 
has a Mittag-Leffier distribution, with Laplace transform E[e~'^^*] = Ei^^—Xt^). Then 
it follows, using (14. 3 p and a simple conditioning argument, that 

poo 

u{t,x)= / E^,[f{Xs);s<rD]Ms)ds 
Jo 

/oo / °o \ 



,n=l 



n=l 

Then, since < E^(-A„t'^) < c/(l + A„t^), we have by ([52D and that 

oo 

||V%||oo < 5^^/3(-A„t^)|(/,^„)|||V^'V^„||oo 

71=1 

OO ,(d+4)/(2a) 

- ^ " l + A^t/^ 



n.=l 

oo 

i+4)/(2a)-l-k 



n=l 



for J = 1,2. Then by f|^ . 

oo 

||V%|U < (cyM)t-''^A(f+^)/(2°)-i- 



n=l 

oo 



n=l 

if A; > (3c/ + 4 - 2a)/(2a). This proves that, when k > -1 + {3d + A) /{2a), u{t,x) is 
in X G K, and hence in D. Consequently, by Remark 13. H the spatial fractional 
derivative A°'/^u{t,x) exists pointwise for x & D, and is a jointly continuous function 
in {t, x). 
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Next we show u{t, x) is CMn t > 0. Let < 7 < 1 A (4/(2a) - 1). By [IHl Equation 

;i7)], 



1 + Xnt^ 

This together with (15.11) and (15. 2p yields that 

5^ |-E,(-A,/) (/,^„)V'n(x)| < J2^Xlt'~'^n'K^^'' 



n=l 



n=l 



< 



n 



{a/d)iy-k+d/{2a)) 



< ct 



7/3-1 



n=l 



Then it follows by a dominated convergence argument that u{t, x) is continuously 
different iable in t > 0, with 



(5.9) 



du{t, x) 



d 



dt 



n=l 



< ct 



7/3-1 



for every x E D. 



Hence by Remark l3. 11 The Caputo fractional derivative d^u{t,x)/dt^ of u{t,x) exists 
pointwise and is jointly continuous in {t,x). Since u(t,x) is a weak solution of (14.81) 
on L^(]R'^; dx), by the above regularity property of u(t,x), it is also a strong solution 
of (Hil). □ 



Remark 5.2. The above proof can be easily modified to show that, if D is a bounded 
regular open subset of M"^ and / G Dom(£|,) for some /c > l+(3(i)/(2a), then ^(t, x) = 
Kx[f{X^J] is a weak solution to the space-time fractional diffusion equation (14. 8p . 
Moreover, the Caputo derivative d^u/dt^ exists pointwise as a jointly continuous 
function in (t, x), and Cdu has a continuous version that equals d^u/dt^ on (0, oo)xD. 



Remark 5.3. The paper [22] solves distributed-order time-fractional diffusion equa- 
tions d^u = An on bounded domains. The distributed-order time-fractional deriva- 
tive is defined by 

where z/ is a positive measure on (0, 1). It may also be possible to extend the results of 
this paper to develop strong solutions and probabilistic solutions for d^u = A°/^u on 
bounded domains. Distributed-order time-fractional diffusion equations can be used 
to model ultraslow diffusion, in which a cloud of particles spreads at a logarithmic 
rate, also called Sinai diffusion [21]. □ 
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Remark 5.4. The fractional Laplacian generates the simplest non-Gaussian stable 
process in M'^. Stable processes are useful in applications because they represent 
universal random walk limits. For random walks with strongly asymmetric jumps, a 
wide variety of alternative limit processes exists, see for example [18]. Because the 
generators of these processes are not self-adjoint, the extension of results in this paper 
to that case remains a challenging open problem. □ 
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